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Abstract 

The possibility of quantum collapse and characteristics of nonlinear localized excitations is ex- 
amined in dense stars with Landau orbital ferromagnetism in the framework of conventional quan- 
tum magnetohydrodynamics (QMHD) model including Bohm force and spin-orbit polarization 
effects. Employ the concepts of effective potential and Sagdeev pseudopotential, it is confirmed 
that the quantum collapse and Landau orbital ferromagnetism concepts are consistent with the 
magnetic field and mass-density range present in some white dwarf stars. Furthermore, the value of 
ferromagnetic-field found in this work is about the same order of magnitude as the values calculated 
earlier. It is revealed that the magnetosonic nonlinear propagations can behave much differently 
in the two distinct non-relativistic and relativistic degeneracy regimes in a ferromagnetic dense 
astrophysical object. Current findings should help to understand the origin of the most important 
mechanisms such as gravitational collapse and the high magnetic field present in many compact 
stars. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. INTRODUCTION 



One of the greatest challenges in astrophysics is to understand the mysterious origin of 
high magnetic field and unique physical processes in many compact stars. There has been 
spectacular breakthroughs in this area of science in the past several years since the pioneering 
works of Chandrasekhar, Bohm, Canuto, Chiu and O'Connell et.al [l-5] on relativistic 
degeneracy and magnetism. We have already observed manu peculiar properties that arise 
due to Pauli exclusion and electron degeneracy in ordinary metallic and semiconductor 
materials such as quantum tunneling, Landau quantization, De Haas- van Alphen effect 
and many others {fJEj]- However, due to known experimental restrictions in astrophysical 
sciences, unfortunately, we are left only with observational and theoretical tools. Initial 
theoretical studies reveal that the relativistic degeneracy js, 9] and magnetism 10-16] can 
lead to quite different equations of state and yet unknown phenomena such as quantum 



collapse 



17, 



18] and Landau orbital ferromagnetism (LOFER) [19l . |20j . A recent comp 



review of different aspects of a dense super-magnetized matter may be found in Refs. 



etc 



21 



22J. Adam [23], using various magnetic models including LOFER, has investigated the 



hydrostatic equilibrium and collapse phenomenon in a white dwarf star and compared the 
Leory with observations. Based on the theoretical calculation of Lee et.al 24j . O'Connell 



25l | has found a stable solution for LOFER sustained self-field which might contribute to 



origin o 
results 



compact astrophysical objects. O'Connell's initial calculation show inconsistent 



26] for the temperatures observed in typical white dwarfs, however, it has been 



shown that the inclusion of electron exchange energy can lead to partial agreement but with 
slightly higher magnetic fields than measured in those objects. 

On the other hand, more recently, Brodin and Marklund have extended the quantum 
magnetohydrodynamics (QMHD) model to include the spin dipole effects and evaluated the 
dynamical behavior of nonlinear waves in an electron-ion plasma. They have reported an 
extreme regimes in which a spin-1/2 quantum plasma exhibits ferromagnetic-type proper- 
ties. In recent years there has been great majority of activities in the field of both quantum 



28 



which reveal 



hydrodynamics (QHD) and quantum magnetohydrodynamics (QMHD) 
the unexpected nonlinear dynamical features in quantum plasmas. Recently, using a classi- 



cal Thomas-Fermi approximation 



40j, 



41] and fully QMHD model 42] it has been confirmed 



that the nonlinear wave properties of a relativistically degenerate plasma fundamentally 
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differs from those of an ordinary degenerate one. Also, other recent studies, employing 



the spin-included QMHD mode 
erties in quantum plasmas 43 



justify the dominant spin effects on nonlinear wave prop- 



511 ] . In this paper, we employ the spin- included quantum 



magnetohydrodynamics model to explore the possibility of Landau orbital ferromagnetism 
and quantum collapse in astrophysical dense objects such as white dwarfs and their pos- 
sible role on nonlinear wave dynamics. The presentation of the paper is as follows. The 
LOFER QMHD model including Bohm potential and the spin-orbit magnetization effects is 
introduced in Sec. [Til and the Localized LOFER density excitations are evaluated and the 
appropriate pseudopotential is found in Sec. IIII1 The numerical results and the comparison 
of findings with the previously reported values are presented in Sec. [TV] Finally, a summary 
is given in Sec. [V] 



II. LOFER QMHD MODEL 



Considering a completely degenerate electron-ion plasma, including the electron spin- 
orbit polarization, quantum tunneling and Lorentz force effects, one may setup the extended 
QMHD equations. For instance, in a center of mass frame of plasma, the continuity equation 
reads as 



Op 
di 



+ v ■ O) = o, 



(i) 



and, the extended momentum equation in cgs units is given as [27] 
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where the field due to the physical current, H, magnetization, M and the sustained field, 
B, are related through the relation; H = B — 47rM. Moreover, the ferromagnetic state of 



plasma is defined by setting H = (i.e. B = 47rM(B)) 52j. Canuto et.al have reported a 



stable LOFER state in strongly magnetized degenerate electron-gas leading to a magnetic 



equation of state of form B = a(p 6 //z e ) 2 / 3 {4, 23, 26], where, p 6 = p/10 6 and p, e is the number 
of nucleons per electron (for instance, p P = 2 for helium atom). The value of proportionality 
constant a is not a priori-evident [23], however, it is known to be sensitively dependent on 
many parameters such as the plasma temperature, electron exchange interactions and the 
maximum magnetic-field oscillation amplitude 20|. It is also noted that the value of a can 
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change due to rapid oscillations in magnetic moment. In what follows, it will be clearly 
evident that our approach and conclusions are independent of the value of this parameter. 
Therefore, the QMHD equations can be closed using a magnetic equation of state such as 
that mentioned above. 

On the other hand, the pressure of a relativistically degenerate magnetized Fermi-gas, 
i.e. P d , (for B/B c < 1, with B c = 4.414 x 10 13 G) is given by [l| 

4 5 

P ~ p d (R) = R (2R 2 - 3) VTTW + Ssinh- 1 ^ , (3) 

where, the Chandrasekhar relativity parameter, R = Pp e /(m e c) = (p/pc) 1 ^ 3 (Pc — 2 x 
10 6 gr /cm 3 ) is a measure of the relativistic degeneracy of electrons and Pp e is the electron 
Fermi relativistic-momentum. It is noted that, in the extreme limits of very small relativity 
parameter, R <C Rch (the value Rch = being the Chandrasekhar critical value of the 
relativity parameter), a dependence Pd oc p 5 / 3 and for corresponding very large values, 
R 3> Rchi a dependence Pd oc p 4//3 is expected for the electron degeneracy pressure. 

For simplicity and without loss of the generality of the problem, we may assume that 
the propagation of magnet oacoustic nonlinear waves takes place along the x direction per- 
pendicular to the high induced magnetic-field B(x,t) along the z-axis. Then, we derive 

dp _,_ dpu _ n 

^ 9X ~ n r , n U) 

dggKg) ,pl , m e c 2 1 dR , ti 2 d r i rl 2 ./7,l 



du I y^du m e c 2 8 
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Note that we have defined a new parameter /3 2 = a 2 {m i jm e c l \ We also define the effective 
quantum potential, the so-called total quantum potential acting on fermions, as 

PR f d R P d {R) JJ3 (3 2 R r— — 
eq = ~ J ~^— dR = "fcT " VTTW - (5) 

The dimensionless set of QMHD equations is thus obtained using the following standard 
scalings 

x — > -^-x, t — > — , p ppo, u — >■ uc s , (6) 



where the normalizing factors, p , u pi = y 'e 2 p /(e m 2 ) and c s = C\fmjm~i denote the equi- 
librium plasma mass-density, characteristic ion plasma frequency and ion quantum sound- 
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speed, respectively. After trivial simplifications, we arrive at the following equations 

dp i dpu _ n. 

at ^~ dx u ' 

§ + up = *i + ^fi^) , (7) 

at. ax her. (it \ */n (1T A I i v f 



dt dx dx dx \ yfp dx 

= - Vi + P 2/3 Rl 

where, the bar notations denote the dimensionless quantities. We have also defined two 
other new parameters, namely, the relativistic degeneracy parameter Rq = (po/Pe) 1//3 and 
the quantum diffraction parameter, H = a/ mi/2m e (hu P i) / (m e c 2 ) interrelated through H = 
eh\/ p c R\jm i 'K I (2m^ 2 c 2 ). The normality notations here and after are avoided for clarity. 
It is realized that the gradient of the effective quantum potential defines the local quantum 
force acting on electrons which is used to describe the quantum-collapse phenomenon and is 
therefore the sum of the quantum-forces due to relativistic electron degeneracy pressure and 



is defined state of 



the relativistic spin-orbit magnetization. The quantum collapse [18 
a plasma for which the effective quantum potential vanishes. 

III. LOFER DENSITY EXCITATIONS 



Since, we are interested in studying the stationary solitary excitations, it is convenient 
to change into the co-moving frame by defining a new variable as £ = x — Vt, and define 
the Eqs. ((7j) in new coordinate. The frame speed, V = u/c s is then the Mach-number of the 
nonlinear magnetoacoustic excitations. Integration of the continuity equation regarding the 
new variable, £, making use of the boundary conditions lim p — 1 and lim u = leads 

£— >±oo £->±oo 

to the relation, u = V (1/p — 1). Using the value for, u, in the momentum equation, with 
p = A 2 , and integrating again with the above-mentioned boundary conditions, reduces the 
Eqs. (JTj) to the following simplified nonlinear differential equation 

= £ (1 - A-r-V> (1 - A - n + ^^-^-^l + ^. ,8) 

Now, multiplying the Eq. flE]) by the quantity dA/dC, and finally integrating with the same 
boundary conditions, gives rise to the well-known energy integral of the form in terms of the 
plasma center of mass density variable 

(dtp) 2 /2 + U(p) = 0. (9) 
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The desired pseudopotential is then given as 



U( P ) 



AnV 2 R 3 - 8nV 2 R 3 p + 3ttR ^1 + R 2 p - 2nR 3 ^/l + R 2 p + R^ 2 p 



-3^0^1 + ^2/3^/3 + 47cV 2 R 3 p 2 + StcR 3 ^TTMp 2 ~ 4R^ 2 p 2 
-QnR 3 p 2 ^Jl + R 2 p 2 / 3 + 3j R4 /3 2 p 7/3 _ ^ pSH^hT 1 Rq + STTpsinh" 1 {R Q p 1 ' 3 ) 



(10) 



The existence of nonlinear localized density structures (solitary waves) requires that the 
following necessary conditions to satisfy, simultaneously 



U(P)\ 



dU(p) 



P =i 



dp 



0. 



d 2 U{p) 



P =i 



dp 2 



< 0. 



(11) 



p=l 



A sufficient condition is to exclude the shock-like structures, i.e., for at least one either 
maximum or minimum there should be at least one nonzero n- value, with U(p m ) = 0, so that 
for all value of p in the range p m > p > 1 (compressive soliton) or p m < p < 1 (rarefactive 
soliton), U (p) is negative. In this way one obtains a potential minimum describing a solitary 
wave propagation given as 

" Pm dp 



± 



-2U(p) 



(12) 



Obviously the two necessary conditions given in Eq. fflT]) is met for the pseudo-potential 
given by Eq. ([101) . The third necessary condition, however, requires the direct evaluation of 
the second derivative of the Sagdeev potential, in Eq. ffTO"]) , at the unstable point, p — 1, 



d 2 U(p) 



dp 2 



P =i 



3H 2 
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2R 
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(13) 



This puts an upper limit on soliton Mach-number for existence of a magnetoacoustic solitary 
excitations the value of which is defined as 



Rq 



RoP 2 

6tt 



(14) 



3^/l + R 2 

The critical Mach-number value obtained above, for = 0, appropriately reduces to the 
value given for the simple unmagnetized Fermi-Dirac plasma [53|]. To this end, we examine 
the sufficient condition, i.e., the existence of the p m values. Thus, we examin the limits of 
the pseudopotential given in Eq. ( FlUj) . i.e. 



lim U (p) 



p-5-0 



> 0, lim U(p) = sgn(/3 2 - 2tt)oo. 



p— >oo 



(15) 
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It is evident that, the solitary waves can only be of rarefactive-type. The possibility of such 



structures has also been reported for relativistic paramagnetic Fermi-Dirac plasma (42] . It 
is also observed that above the critical value of (3 cr = y/2n nonlinear excitations (periodic or 
solitary) are forbidden. In the next section we will again encounter this critical value and 
its prominent definition in the context of quantum collapse in orbital-ferromagnetic plasma. 

IV. DISCUSSION 

Figure 1 depicts the thick curve where the effective quantum potential vanishes, i.e., 
ty eq = 0, and the filled region, where, the magnetosonic solitary excitations are stable. 
It is evident from Eq. (fl5|) that no solitary or periodic excitations exist above a critical 
value of (3 cr > y/2ir. It is also noted that the quantum collapse in LOFER plasma state 
resides in the region where the nonlinear excitations disappear. The quantum collapse 



occurs for j3 = ^/2ty/R(1 + R 2 ) 1 ^, or equivalently, for the induced field strength of B = 
C\/27rm e /mj_R 3//2 (l + i? 2 ) 1 / 4 . It is noted that the (3 cr = \[2tx value corresponds to the value 
of a ~ 1.754 x 10 9 , which is slightly higher than the value reported by O'Connell and Roussel 



25l | but lower than the value reported by Pohl and Schmidt-Burgk 20j . The vertical line in 
Fig. 1 corresponds to Chandrasekhar critical relativity parameter value of R = a/2 which 
corresponds to the magnetic field strength of B ~ 3.882 x 10 9 G in the collapsing LOFER. 
The corresponding magnetic field strength can be found in a highly magnetized white dwarf. 
For a typical inertial confined dense plasma with electron density of n e ~ 10 23 cm -3 it is 
easily found that, in laboratory LOFER, the quantum collapse can occur for the magnetic 
field of B = 9.737 x 10 6 G which might be realized in future. As it was mentioned earlier, 
parameter a and consequently /3 contains a rapidly oscillating part changing as the successive 
Landau magnetic- levels are filled/unfiiled by the electrons, so that there is a possibility 
for the LOFER state to oscillate back and forth between nonlinear excitation region and 
quantum collapse region, when the /3 is around the critical-value. This effect may lead to 
the magnetosonic pulsation when the density is very high or when the the Landau levels are 
very close to each other. 

In Fig. 2 we show the volume in V-Rq-(3 space where the localized nonlinear excitations 
exist. It is clearly remarked that, for a fixed value of (3, there always exists a lower cut-off 
degeneracy parameter value Ro m below which the solitary excitations are absent. It is also 



7 



observed that, as the value of (5 increases lower Mach-numbers will be available to localized 
nonlinear excitations and conversely for the fixed values of /5 as the degeneracy parameter 
i?o increases the higher Mach-numbers will be available to solitary excitations. On the other 
hand, the pseudopotential potential and their shape variations (when only one fractional 
plasma parameter is changed and others are fixed) are shown in Fig. 3. It is revealed that, 
the change in the pseudopotential profile with change in plasma parameters is more sig- 
nificant for non-relativistic degeneracy cases (left column) rather than the relativistic ones 
(right column). This feature is further revealed in Fig. 4, where the rarefactive localized 
excitation are depicted. From Fig 4. (which borrows the same parameters from Fig. 3) 
is easily concluded that there is fundamental differences between non-relativistic and rela- 
tivistic localized nonlinear excitations. Fig. 4(a) shows that for the non-relativistic LOFER 
plasma the soliton amplitude/width increases/decreases as the degeneracy parameter is in- 
creased for fixed value of f3, while, Fig. 4(b) reveals that in the relativistic case the amplitude 
is almost unaffected by this variation and only the soliton with increases with increase in 
the value of relativistic degeneracy. Concerning the variation of (3 (or equivalently the value 
of induced magnetic field) one observes from Fig. 4(c) that the increase in the value of (3 
for a fixed value of the degeneracy parameter, R , opposite to the case of Fig. 4(a) causes 
a decrease/increases in the soliton amplitude/width. Finally, Fig. 4(d) reveals that the 
change in the parameter /3 for a fixed value of degeneracy parameter value in the case of 
relativistic degeneracy leaves the soliton amplitude and width almost invariant opposite to 
the case of non-relativistic degeneracy. The mentioned fundamental differences between the 
cases of non-relativistic and relativistical LOFER plasma degeneracy mark distinct features 
similar to the ones presented in Ref. 42] and can help understand the underlying physical 
nonlinear processes in stellar dense plasmas. 



V. SUMMARY 



We investigated the possibility of a quantum collapse and nonlinear localized excita- 
tions in a Landau orbital ferromagnetic dense stellar plasma in the framework of quantum 
magnetohydrodynamic model using the concepts of effective potential and Sagdeev pseu- 
dopotential. Our findings confirm that a quantum collapse is possible for magnetic field 
strengths present in some white dwarf stars and the value of induced field found here are 
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consistent with the ones reported earlier. Our investigation also reveals that the magne- 
tosonic nonlinear propagations behave rather differently in the two distinct non-relativistic 
and relativistic degeneracy regimes. Current findings can help better understand the origin 
of the most important mechanisms such as collapse and high magnetic field observed in some 
compact astrophysical objects. 
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FIGURE CAPTIONS 



Figure- 1 

The zero effective potential (quantum collapse) curve (thick-curve) and the region (filled- 
region), where the magnetosonic localized excitations occur for Landau orbital ferromagnetic 
(LOFER) plasma. The vertical line denotes the Chandrasekhar critical relativity parameter 
value (R=y/2). 

Figure-2 

Figure 3 shows a volume in V-(3-Ro space in which a solitary excitation exists. 
Figure-3 

(Color online) The variations of pseudopotential depth and width for rarefactive solitary 
nonlinear magnetosonic waves with respect to change in each of three independent plasma 
fractional parameter, normalized soliton-speed, V, the ferromagnetic parameter, (3, and the 
relativistic degeneracy parameter, R , while the other three parameters are fixed. The 
dash-size of the potential curves increase with regard to increase in the varied parameter. 

Figure-4 

(Color online) The variations of localized magnetoacoustic density-excitation profiles with 
respect to change in each of three independent plasma fractional parameter, normalized 
soliton-speed, V, the ferromagnetic parameter, j3, and the relativistic degeneracy parameter, 
R , while the other three parameters are fixed. The thickness of the profile increase with 
regard to increase in the varied parameter. 
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Quantum Collapse (i^ e[1 = 0) and Soliton Stability Region 
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Non-Relativistic Degeneracy Relativistic Degeneracy 
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